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Abstract. We present the results of a photometric multisite campaign on the 6 Scuti Pre-Main-Sequence star IP Per. Nine 
telescopes have been involved in the observations, with a total of about 190 hours of observations over 38 nights. Present data 
confirms the multiperiodic nature of this star and leads to the identification of at least nine pulsational frequencies. Comparison 
with the predictions of linear non-adiabatic radial pulsation models allowed us to identify only five of the nine observed frequen¬ 
cies, and to constrain the position of IP Per in the HR diagram. The latter is in good agreement with the empirical determination 
of the stellar parameters obtained by [Miroshnichenko et al. (2001)] An initial interpretation of the observed frequencies using 
the Aarhus non-radial pulsation code suggests that three frequencies could be associated with non-radial (1=2) modes. 

Finally, we present new evolutionary and pulsation models at lower metallicity (Z=0.008) to take into account the possibility 
that IP Per is metal deficient, as indicated by | Miroshnichenko eFalT (200 T) | 

Key words, stars: variables: 6 Scuti - stars: oscillations - stars: pre-main sequence - stars: fundamental parameters - individual 
IP Per 


1. Introduction 

Intermediate mass (M >1.5 Mq) Pre-Main-Sequence (PMS) 
stars are known as Herbig Ae/Be stars l |Herbig 1960^ . This 
class of stars is characterized by spectral type A or B 
with emission lines, an infrared excess due to hot or cool 
circumstellar dust or both, and luminosity class III to V 
dWaters & Waelkens I998I I. Herbig Ae/Be are also well known 
for their photometric and spectroscopic variability on time 
scales of minutes to years mainly due to photospheric activ¬ 
ity and interaction with the circumstellar environment (see e.g. 
ICatala 20031 1. However, the fact that these young stars dur¬ 
ing their contraction towards the Main-Sequence (MS) move 
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across the pulsation instability region of more evolved stars 
has prompted the suggestion that at least part of the activ¬ 
ity could be due to stellar pulsation Isee lBaade & Stahl 19891 
|Kurtz & Marang 1995t . 

The possible presence of pulsators among Herbig Ae/Be 
stars is particularly attractive since the precise observables 
which can be measured, i.e. the pulsation frequencies can, in 
principle, allow us to test evolutionary models by constraining 
the internal structure using asteroseismological techniques. 

The existence of pulsating Herbig stars was originally sug¬ 
gested by |Breger (1972)| who discovered two candidates in the 
young open cluster NGC 2264. This initial finding was con¬ 
firmed by subsequent observations of 5 Scuti-like pulsations 
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Table 1. Journal of the observations. Note that the data coming 
from Sierra Nevada and San Pedro Martir (SPM) have been 
transformed in the Johnson system as discussed in section 2.2. 


HJD-2450000 
start (days) 

HJD-2450000 
end (days) 

Duration 

(hours) 

Filter 

2545.585 

Loiano CCD (Italy) 

2545.659 

1.8 

B 

2547.599 

2547.681 

2.0 

B 

2548.547 

2548.611 

1.5 

B 

2549.503 

2549.623 

2.9 

B 

2628.347 

2628.627 

6.7 

B 

2629.286 

2629.453 

4.0 

B 

2630.454 

2630.513 

1.4 

B 

2654.252 

2654.332 

1.9 

B 

2976.875 

Fairborn-APT (USA) 

2976.892 

0.4 

BV 

2977.851 

2977.936 

2.0 

BV 

2978.721 

2978.933 

5.1 

BV 

2983.573 

2983.920 

8.3 

B 

2984.711 

2984.917 

4.9 

BV 

2987.765 

2987.902 

3.3 

BV 

2977.170 

SOAO (Korea) 

2977.333 

3.9 

BV 

2980.910 

2981.063 

3.7 

V 

2983.007 

2983.191 

4.4 

BV 

2988.267 

OSN (Spain) 

2988.500 

5.6 

uvby 

3047.364 

3047.472 

2.6 

uvby 

3048.299 

3048.459 

3.8 

uvby 

2946.330 

Loiano TTCP (Italy) 

2946.703 

9.0 

BV 

2947.287 

2947.441 

3.7 

BV 

2973.393 

2973.636 

5.8 

BV 

2970.631 

SARA (USA) 

2970.756 

3.0 

B 

2970.494 

SPM (Mexico) 

2970.818 

7.8 

uvby 

2972.495 

2972.841 

8.3 

uvby 

2973.493 

2973.775 

6.8 

uvby 

2974.487 

2974.585 

2.4 

uvby 

2977.485 

2977.847 

8.7 

uvby 

2972.974 

BAG (China) 

2973.383 

9.8 

V 

2974.032 

2974.146 

2.7 

V 

2974.939 

2975.382 

10.6 

V 

2976.940 

2977.078 

3.3 

V 

2984.196 

2984.292 

2.3 

V 

2977.506 

Teide-OGS (Spain) 

2977.603 

2.3 

V 

2966.462 

Serra la Nave (Italy) 

2966.62 

3.9 

BV 

2967.337 

2967.61 

6.7 

BV 

2973.348 

2973.61 

6.4 

BV 


in the Herbig Ae stars HR5999 (|Kurtz & Marang 1995J and 
HD104237 JDonati et al. 19^ . 

This empirical evidence stimulated the hrst theoretical in¬ 
vestigation of the PMS instability strip based on non-linear 
convective hydrodynamical models (Marconi & Palla 1998). 
As a result, the topology of the PMS instability strip for the hrst 
three radial modes was identihed. IMarconi & Palla (1998)| also 
pointed out that the interior structure of PMS stars entering the 
instability strip differs signihcantly from that of more evolved 
Main Sequence stars (with the same mass and temperature), 
even though the envelopes structures are similar. This property 
was subsequently conhrmed by |Suran et al. (2001)| who made 
a comparative study of the seismology of a 1.8 Mq PMS and 
post-MS star. |Suran et al. (200 1)| found that the unstable fre¬ 
quency range is roughly the same for PMS and post-MS stars, 
but that some non-radial modes are very sensitive to the deep 
internal structure of the star. In particular, it is possible to dis¬ 
criminate between the PMS and post-MS stage using differ¬ 
ences in the oscillation frequency distribution in the low fre¬ 
quency range (g modes, see also |Templeton & Basu 2003) . 

Up to now new observational programs have been car¬ 
ried out by various groups. The current number of known 
or suspected candidates amounts to about 30 stars (see the 
updated list at http://ams.astro.univie.ac.at/pmsxorot.php and 
the reviews by IZwintz et al. 20041 IMarconi & Palla 20041 and 
IMarconi et al. 20041 . However, only a few stars have been stud¬ 
ied in detail, so that the overall properties of this class of vari¬ 
ables are still poorly determined. 

In this context, our group has started a systematic mon¬ 
itoring program (see IMarconi et al. 20()lT jRipem^taL^^0^ 


IPinheiro et al. 200^ [Ripepi et al. 2003| IB ernaberetar2004 


of Herbig Ae stars with spectral types from A to F2-3 with 
the following aims; 1) to identify the largest number of 
pulsating objects in order to observationally determine the 
boundaries of the instability strip for PMS d Scuti pulsation; 
2) to study in detail through multisite campaigns selected 
objects showing multiperiodicity (see IMarconi et al. 20011 
RiDeDi^etaL^2002 jPinherro et al. 20031 [Ripepi et al. 2003| 


Bernabei et al. 20041 . The multiperiodic pulsators are potential 


candidates for future asteroseismological analysis. 

During this project our attention turned to the star IP Per, al¬ 
ready listed as Herbig Ae star by [The et al. (1994) with spectral 
type A3e, and studied in detail bylMiroshmchento et al. (2001) 


both photometrically and spectroscopically. The main proper¬ 
ties of this interesting object are: 


- IP Per is a typical UX Ori type star showing photometric 
variations with an amplitude of ~0.3 mag and with duration 
of the minima of 10-50 days. 

- the fundamental stellar parameters are: spectral type A7, 

Teff ~ 8000 + 100 K, logg ~ 4.4 + 0.1, = 2.0 + 1.0 

km s“^ and [M/H] ^ -0.41 +0.1. Thus, IP Per is a dwarf 
metal-poor star. 

- its radial velocity and proper motion suggest that it most 
likely belongs to the Per OB2 association (D ^ 300 pc, 
Ide Zeeuw et al. 19991 . At Z) ^ 300 pc and logL/L q-1.0, 
IP Per falls onto the ZAMS at a mass of 1.8 Mq. 
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Table 2. List of telescopes and instruments involved in the mul¬ 
tisite campaign. SCP means Single Channel Photometer; TCP 
means Three Channel Photometer. To give an idea of the rela¬ 
tive photometric precision of each site, the last column shows 
the white noise in the Fourier space (amplitude) for each site. 
All the calculations are in the V filter except for Loiano (2002) 
and SARA sites for which we have only B filter data. 


Observatory 

Telescope 

Instrument 

noise 

(mmag) 

Loiano (Italy) 

1.5m 

BFOSC- 2002 

0.2 

Loiano (Italy) 

1.5m 

TTCP - 2003 

0.2 

BAO (China) 

0.85m 

TCP 

0.4 

SPM (Mexico) 

1.5m 

uvby Phot. 

0.6 

SARA (USA) 

0.9m 

CCD 

1.5 

Teide (Spain) 

1.0m OGS 

CCD 

0.6 

Fairborn (USA) 

0.75m T6 

SCP 

0.8 

OSN (Spain) 

0.9m 

uvby Phot. 

0.6 

Serra la Nave (Italy) 

0.9m 

SCP 

1.3 

SOAO (Korea) 

0.6m 

CCD 

0.7 


In the following, we present the results of a photomet¬ 
ric multisite campaign to study in detail the pulsational prop¬ 
erties of IP Per. In Section 2. we describe the observations 
and data reduction techniques; in Section 3. we discuss the 
frequency analysis; Section 4 presents a comparison with 
theoretical models. A brief discussion concludes the paper. 
Preliminary results of our observations have been presented in 
[Ripepi & Marconi (2004)| and |Ripepi et al. (2004)| 
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Fig. 1. Light curves for the 2002 B dataset (Loiano). Note that 
SB = Bvar - Bcomp- The solid line represents the fit to the data 
including the six significant frequencies found for this dataset 
(see Sect.3). 


duced using standard IRAF routines. Aperture photometry was 
carried out by using the DAOPHOTII package (IStetson 19871 . 

The multisite campaign in 2003 was carried out by using 
a variety of instruments, as described in some details in the 
Appendix. 


2. Observations and data reduction 

The campaign on IP Per was conducted in two parts; 1) single 
site observations during winter 2002/2003 (Loiano telescope); 
2) multisite campaign during winter 2003/2004, involving 9 
different telescopes/instruments, as described in Table |2] and 
more in detail in the next section. A total of about 190 hours of 
observations spanning 38 nights have been gathered. A detailed 
log of the observations is shown in Tabled 

As for the comparison star, we typically used the star TYC 
2359-802-1 (03* 41'” 15.36" +32° 23' 49.7" (2000) B7-=12.029 
mag, V 7 -=l 1.299 mag). It has been checked on the basis of 
CCD data obtained in Loiano in 2002, see below). HD 278941 
(03* 41”' 01.67" +32° 07' 06.l" (2000) Bt^9.957, Vr=9.576 
mag, A5) was the check star. 

2.1. Details on single site observations 

The exploratory run on IP Per was carried out in 2002/2003 
with the Loiano L5m telescope (Italy). The observations has 
been gathered in the B filter by using the BFOSC instrument 
(http;//www.bo.astro.it/loiano/observe.htm#manuals i equipped 
with a CCD EEV 1300x1340 pixel (R.O.N.=L73 e-/pixel; 
GAIN=2.13 e-/ADU). The pixel scale was 0.58 arcsec/pixel, 
for a total field of view of 12.6 x 13 acrmin^. Sky flats, dark 
and bias exposures were taken every night. All data were re- 


2.2. Merging the data from different sites 

In order to use simultaneously all the data gathered during 
the campaign we had to convert the data coming from Sierra 
Nevada and San Pedro Martir sites from Stromgren uvby Alters 
to the Johnson B and V ones. This was straightforward for the 
V because bV ^ by (where bV means Vyarmbk - Vcomparison), 
whereas to obtain bB we used the linear approximation by 
[Warren & Hesser (1977)| (equation 5), obtaining bB by + 
1.668(5(17-y). 

The different datasets from the various telescopes were in¬ 
spected carefully, and scattered points due to bad weather (very 
frequent for data obtained with single channel photometers 
data) or other causes have been removed. 

Due to the differences of the various instruments and Alters 
used, and in order to prevent problems with zero point differ¬ 
ences between different datasets, as well as to prepare the data 
for Eourier analysis, we decided to detrend the data to a com¬ 
mon average zero value. 

In total, during our study of IP Per we obtained three dif¬ 
ferent time series: 

1. B Alter data from Loiano 2002/2003 observations (2002 B 
dataset here after) 

2. B Alter data from multisite campaign 2003/2004 (2003 B 
dataset here after) 
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Fig. 2. Left panels; light curves for the 2003 V dataset. Note that 6V - Vvar - Vcomp- Right panels; light curves for the 2003 B 
dataset. Note that dB = Bvar - Bqomp 


3. V hlter data from multisite campaign 2003/2004 (2003 V 

dataset here after) 

The light curve for the 2002 B dataset is shown in Fig. [J 
while the left and right panels of Fig.|3show the light curves for 
the 2003 V and B datasets respectively. In order to appreciate 
the quality of the V data taken in 2003 we report in Fig. 13 a 
zoom of Fig. 13 (left panel) including only the Loiano TTCP 
data (V filter), which is representative of our highest quality 
photometry. 


2.3. uvby Hr measurements for IP Per 


During our uvby measurements collected at Sierra Nevada 
Observatory, two brighter comparison stars C1 =HD 22418 and 
C2=HD 21913 were used for purposes of calibration. In ad¬ 
dition, a few Hp observations were also obtained. Then, in¬ 
strumental magnitude differences were obtained relative to 
Cl. To transform these instrumental differences into the stan¬ 
dard uvbyp system, we have followed the procedure described 
in |Rodrfguez et al. (1997)| Thus, the transformation equations, 
obtained using a set of 19 standards stars selected from the list 
of [ Crawford & Mander (196^ and fCrawford & Barnes (1970) 
were used. 

Next, the absolute standard uvbyp indices of IP Per, Cl 
and C2 were obtained following the method described in 
[Rodriguez et al. (2003 )| using Cl and C2 as zero-points. The 
results are listed in Table|3together with those listed in the bib¬ 
liography for the comparison stars. The error bars in this table 
represent the standard deviations of magnitude differences rela¬ 
tive to Cl. As seen, our results are in very good agreement with 
the values found in the homogeneous catalogue of Olsen (1996, 


Table 3. uvbyp indices obtained for IP Per and comparison 
stars. The pairs below the star names are the number of points 
collected for each object in uvby and /?, respectively. The num¬ 
bers below the measures show the errors (in mmag). The values 
given by Olsen (1996, private communication) are listed in the 
bottom part. 


Star 

V 

b-y 

mi 

Cl 

P 


(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

IP Per 

10.374 

0.239 

0.144 

0.834 

2.763 

(49,5) 

6 

6 

6 

15 

9 

C1=HD22418 

6.963 

0.280 

0.163 

0.420 

2.646 

(49,6) 

4 

3 

4 

9 

4 

C2=HD21913 

7.627 

0.297 

0.136 

0.473 

2.640 

(21,5) 

4 

3 

4 

9 

5 

C1=HD22418 

6.962 

0.281 

0.167 

0.419 

2.645 

C2=HD21913 

7.628 

0.296 

0.132 

0.473 

2.640 

private communication). Similar results can be found in the 

bylOlsen E.H. 

(l983)|and 

Hauck B. & Mermilliod (1998) 


3. Frequency Analysis 

The frequency analysis was performed using the period04 
package ( jLenz & Breger 2005^ , based on the Fourier transform 
method. For a better interpretation of the results, we have first 
calculated the spectral window (SW) for each dataset. The re¬ 
sult is shown in Fig.0]where from top to bottom we report the 
SW for the three time series identified in the previous section 
(see labels in the figure). The SW was used as a diagnostic to 
distinguish between real and spurious frequencies. 
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Table 4. Frequencies, amplitudes and confidence levels for the three datasets analysed in this paper. The error bars on the 
frequencies are +0.11 c/d, ±0.15 c/d and ±0.25 c/d for the datasets 2003 V, 2003 B and 2002 B datasets, respectively. Note that 
for the 2003 B and 2002 B datasets the frequencies are not listed by decreasing amplitude, but by correspondence with the 2003 
V dataset frequencies (see text). 


2003 V dataset 2003 B dataset 2002 B dataset 



Frequency 

(c/d) 

Amplitude 

(mmag) 

confidence 

(%) 

Frequency 

(c/d) 

Amplitude 

(mmag) 

confidence 

(%) 

Frequency 

(c/d) 

Amplitude 

(mmag) 

confidence 

(%) 

/l 

22.89 

1.9 

99.9 

22.89 

3.1 

99.9 

22.88 

3.3 

99.9 

/2 

34.60 

1.5 

99.9 

34.82 

2.1 

99.9 

34.64 

3.3 

99.9 

h 

30.45 

1.8 

99.9 

30.45 

1.8 

99.9 

30.48 

3.2 

99.9 

h 

48.23 

1.6 

99.9 

48.45 

2.0 

99.9 

48.23 

1.9 

99.9 

h 

28.79 

1.5 

99.9 




27.73 

1.9 

99.9 

h 

23.99 

1.3 

99.9 







fl 

9.30 

1.3 

99.9 







h 

41.11 

1.2 

99.9 

42.08 

1.9 

99.0 

42.27 

2.3 

99.9 

h 

52.00 

1.1 

99.9 









HJD-a453000 

Fig. 3. Data collected in Loiano with the TTCP (V filter) dur¬ 
ing Nov./Dec. 2003 (dots). The solid line show the fit to the 
data with 9 frequencies obtained from the period analysis (see 
Sect.3). 


Each data set described in the previous Section has been 
analysed separately. Figures |3 E] 0 show the Fourier trans¬ 
form for the datasets 2002 B, 2003 B and 2003 V respec¬ 
tively. Here, in each panel the peak with largest amplitude is 
selected and then removed, obtaining a new spectrum shown 
in the following panel. The last panel shows the periodogram 
after the prewithening with all the significant frequencies. The 
solid, dashed and dotted lines show the 99.9%, 99% and 90% 
confidence levels calculated following the widely used recipe 
by |Breger et al (199^ and |Kuschnig et al. (1997)| The eiTor on 
the measured frequencies (apart from the ± 1 c/d alias) can 
be roughly estimated from the FWHM of the main lobe in 
the spectral window (see lAlvarez et al. 199^ and references 


therein). As a result we found A/ ~ 0.25 c/d, 0.15 c/d and 0.11 
c/d for the datasets 2002 B, 2003 B and 2003 V respectively. 

The frequencies found for the three datasets are summa¬ 
rized in Table 0 In order to discuss in detail the results of the 
frequency analysis, let us take as reference the frequencies ob¬ 
tained with the 2003 V dataset which is the best one. Then we 
find; 

- /i = 22.89 c/d, /2 = 34.60 c/d, /a = 30.45 c/d, = 48.23 
c/d are present, within the quoted errors, in every dataset. 

- f 5 - 28.79 c/d is the +1 c/d alias of the / = 27.73 c/d 
frequency found in the 2002 B dataset. Due to the better 
duty cycle of the V data we shall retain - 28.79 c/d for 
comparison with the models. 

- h - 41.11 c/d is probably the -1 c/d alias of / = 42.08 
c/d in the 2003 B dataset and / = 42.27 c/d in the 200V B 
dataset. As before, we shall use the frequency obtained in 
the V filter. 

- /6=23.99 c/d, fi-9.30 c/d and fg - 52.00 c/d are found 
only in the V filter probably because of the better S/N of 
this dataset. 

In summary, the multisite campaign allowed us to identify 
up to 9 significant frequencies for IP Per. The fit to the observed 
data points obtained with period04 including the nine is shown 
in Fig0 Five of these are present in each dataset. In particular 
we confirm all the six frequencies (within the errors and taking 
into account the 1 c/d alias) found during the single site CCD 
study (i.e. the 2002 B dataset) of IP Per. 

4. Comparison with theoreticai modeis 

The results presented in the previous Section can be used 
to constrain the intrinsic stellar properties of IP Per and in 
particular its mass and position in the HR diagram, through 
comparison with stellar pulsation models. Using a linear 
non-adiabatic pulsation code (see Marconi & Palla 1998, 
IMarconi et al. 20041 for details) we could not reproduce all 
the observed frequencies. In fact, we can recover at most 
5 of the 9 observed frequencies for M - 1.77 ± O.OIMq, 
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Fig. 4. Spectral window in amplitude for the three data sets 
analysed in this paper (see labels in the figure). Note in par¬ 
ticular that the 1 c/d alias is reduced to about 50% for the 2003 
V dataset. 


Fig. 6. As in Fig.|3but for the 2003 B dataset. 
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Fig. 5. Frequency analysis for the 2002 B dataset. The solid, 
dashed and dotted lines show respectively the 99.9%, 99% and 
90% significant levels. In each panel, one peak (i.e. the la¬ 
belled frequency) is selected and removed from the time se¬ 
ries and a new spectrum is obtained. The last panel shows the 
periodogram after the prewithening with all the significant fre¬ 
quencies. 
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Fig. 7. As in Fig.|3but for the 2003 V dataset. 


logL/Lo = 0.992 + 0.003, logTes = 3.887 + 0.002. This so¬ 
lution corresponds to a radial pulsation model which simul¬ 
taneously oscillates in the first (/i), second (/s), third (/ 2 ), 
fifth (/ 4 ) and sixth (fg) overtones. Its position in the HR di¬ 
agram is shown in Fig. together with the predicted insta¬ 
bility strip by [Marconi & Palla (1998)| and the PMS evolu¬ 
tionary tracks computed for the labelled stellar masses with 
the FRANEC stellar evolution code ( IChieffi & Straniero l989j 
ICastellani et al. 1999> . The I.IIMq PMS track is represented 
by the dotted line. 
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Effective Temperature Log (K) 


Fig. 8. Position in the HR diagram of the best-fit radial pul¬ 
sation model. The shaded region is the predicted instabil¬ 
ity strip by lMarconi & Palla (1998)| with solar chemical com¬ 
position (Z-0.02, Y=0.28). Solid lines are the PMS evolu¬ 
tionary tracks computed for the labelled stellar masses with 
the FRANEC stellar evolution code (Chieffi & Straniero 19891 
ICastellani et al, 19991 . The dashed lines represent the I.IIMq 
PMS track with solar chemical composition. 


In order to reproduce the observed frequencies that 
were not found from the radial analysis and, at the same 
time to investigate the possibility, strongly supported by 
empirical evidence on other PMS 6 Scuti stars (see e.g. 
IBalona et al. 2002V that non-radial modes are also present in 
IP Per, we have attempted an asteroseismological interpreta¬ 
tion of the data using the Aarhus adiabatic non-radial pulsa¬ 
tion code (http;//astro.phys.au.dk/~jcd/adipack.n/). A prelimi¬ 
nary application of this code to the PMS evolutionary structure 
corresponding to the best-fit radial pulsating model seems to 
suggest that /e and /s are associated to non-radial modes 
with / = 2. As for fi, neither the radial nor the non-radial anal¬ 
yses are able to match the observed value for the selected stellar 
parameters. 

4.1. Comparison with post-MS puisation modeis 

Both spectroscopic measurements and PMS model predictions 
seem to suggest that the position of IP Per in the HR diagram is 
near the MS, in a region where PMS and post-MS evolutionary 
tracks, at fixed stellar mass, are known to intersect and remain 
quite close to each other. In order to investigate the effect of the 
assumed evolutionary status on the predicted frequencies, we 
have taken into account the post-MS evolution of stellar mod¬ 
els with masses ranging from 1.7 to 1.8 Mq, computed with 
the FRANEC code. As a result, we find that the post-MS evo¬ 
lutionary model located at the same luminosity and effective 
temperature of our PMS best fit model (see Fig.|^, has a stellar 


mass of 1.73 Mq, i.e. only slightly lower than the PMS counter¬ 
part (1.77 Mq). As a consequence the structure of this post-MS 
model produces periods slightly longer than the best fit PMS 
solution, but differences are of the order of minutes for both 
radial and non-radial p modes. The 1.73 Mq post-MS model, 
reproducing the same number of frequencies as the PMS one, 
would require a slightly different luminosity and effective tem¬ 
perature, namely logL/Lo = 0.988, logTeff = 3.888. The val¬ 
ues of the predicted p mode periodicities are in this case very 
similar to the PMS ones. However, we know that the small 
differences in the non-radial frequencies can produce evident 
changes in the large and small frequency separations, in partic¬ 
ular for / = 2, as extensively discussed by |Suran et al. (200 1)| 
This occurrence would in principle allow to verify the PMS na¬ 
ture of IP Per if very accurate and long time based observations 
(e.g. from space) were available. 


4.2. Effect of the metai abundance 


All the above evolutionary and pulsational analysis is 
based on models with solar chemical composition, 
namely Z=0.02, Y=0.28. If the evidence pointed out by 
Miroshnichenko et al. (2001) that IP Per has significantly 


lower metal abundance (Z ^ 0.008) is confirmed, the above 
theoretical interpretation would have to be modified, as 
discussed in the following. 


4.2.1. The instability strip 

We computed new nonlinear radial pulsation models, with the 
same code and the same numerical and physical assumptions 
as in Marconi & Palla (1998), but with Z=0.008, Y=0.25. The 
resulting instability strip for the first three radial modes is re¬ 
ported in Figl^ (dashed lines). As shown in this figure, where 
the instability strip for Z=0.02 is also plotted for comparison, 
the metallicity effect on the theoretical boundaries is rather 
small (200 K at most for each luminosity level) and decreases 
toward the higher luminosities, with a global shift of the in¬ 
stability region toward higher effective temperatures, as Z de¬ 
creases from Z=0.02 to Z-0.008. 


4.2.2. Pulsational frequencies 

We also computed new PMS evolutionary models at Z=0.008, 
Y=0.25 with the FRANEC code and the same assumptions 
used for solar chemical composition. Then we computed lin¬ 
ear nonadiabatic radial pulsation models* along the Z=0.008 
PMS tracks, in order to match the observed frequencies of IP 
Per. The resulting best fit radial model for Z=0.008 Y=0.25, 
again pulsating in five radial modes with /i, / 2 , ft, fs and fg, 
corresponds to a PMS structure with M - 1.54 + O.OIMq, 
logL/Lg = 0.991 +0.003, log Tefr = 3.896 + 0.002. Its position 
in the HR diagram is shown in Fig. [^together with the pre- 

’ We did not compute nonradial models because our version of the 
Aarhus nonradial code assumes solar metallicity. The extension of 
nonradial computations to other chemical compositions will be ad¬ 
dressed in a future work. 
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Effective Temperature Log (K) Effective Temperature Log T^ (K) 


Fig. 9. Position of the instability strip for the three radial modes 
at Z=0.008 (dashed lines) compared with the one at solar 
metallicity (solid lines). The long dashed line represent the 
birthline by |Palla & Stabler (1993)| 

dieted instability strip for Z=0.008 and the corresponding PMS 
evolutionary track. Inspection of this hgure clearly indicates 
that the predicted position in the HR diagram of IP Per, assum¬ 
ing Z=0.008, is in good agreement with the empirical spec¬ 
troscopic determination by |Miroshnichenko et al. (20Ql)j (hlled 
circle). 

However, the observationally derived metal deheit reflects 
only the current state of the atmosphere and may be due to re¬ 
cent effects, such as interaction with the circumstellar medium 
(see [Gray & Corbally 1998^ . In this case it would represent 
only a surface effect. As the PMS structure of our best ht model 
is characterized by a radiative envelope, such a contamination 
should not have any effect on the pulsation properties which 
involve deeper layers across the Hydrogen and Helium ioniza¬ 
tion regions.On the other hand the different abundances in the 
surface layer would affect the position in the HR diagram due 
to the opacity variation. A detailed modeling of this case is 
beyond the purposes of the present paper. At the same time a 
more accurate determination of the metal poor nature of IP Per 
would be important. 


5. Conclusions 

A total of about 190 hours of observations obtained during 38 
nights at 9 different telescopes on the PMS 6 Scuti star IP Per 
have been presented. The Fourier analysis of this data set con- 
hrms the multiperiodic nature of this pulsator: we have identi- 
hed nine frequencies of pulsation which are signiheant on the 
basis of the |Breger et al (1993)| and |Kuschnig et al. (1997)| cri- 
teria. 

Comparison with the predictions of linear non-adiabatic ra¬ 
dial pulsation models allows us to recover only five out of the 


Fig. 10. Best ht radial pulsation model for Z=0.008 in the 
HR diagram (hlled triangle) as compared with the spectro¬ 
scopic measurement by Miroshnichenko et al. (2001) The cor¬ 
responding PMS evolutionary track and instability strip are 
overplotted. The track for 1.77 Mq, Z=0.008 is shown for com¬ 
parison with Fig.|^ 


nine observed frequencies for a 1.77 Mq model. Non-radial 
pulsation is also present in this star. A preliminary interpre¬ 
tation of the observed frequencies through the Aarhus non- 
radial code, applied to the evolutionary structure of the 1.77 
Mq model reproducing /i, / 2 , fs and fg with radial modes, 
seems to indicate that f^, /e and f% are associated with non- 
radial modes with 1-2. Specihe post-MS evolutionary and 
pulsational models were computed in order to investigate the 
dependence of radial and non-radial output frequencies on the 
assumed evolutionary status. The resulting post-MS solution 
has similar stellar parameters and p mode frequencies. 


Finally the possible effect of the metal poor nature of IP Per 
detected by Miroshnichenko et al. (2001) on both pulsation and 
evolutionary properties is discussed. We hnd that if the metal¬ 
licity of IP Per is as low as Z=0.008 the best fit radial model 
has a significantly lower mass than the case at solar chemical 
composition but the pulsation characteristics are similar.Also 
the estimated position of IP Per in the HR diagram appears to 
be in good agreement with the independent determination by 
Miroshnichenko et al. (2001) Whether the low metallicity is a 
property only of the surface layers or represents a systematic 
deheit throughout the interior, as we have assumed in our mod¬ 
eling, should be clarihed before hnal conclusions on the stellar 
parameters of IP Per can be reached. 
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Appendix A: Details on the instrumentation used 

In this appendix we describe in some detail the combinations 
telescope/instrument used during the 2003 multisite campaign 
and summarized in Tab.|2] An indication on data reduction pro¬ 
cedures is reported too. 

1. The Loiano (2003) observations were obtained 

with the 1.5m telescope equipped with the 
Toppo Three Channel Photometer (TTCP, see 

http://www.na.astro.it/%)7Esilvotti/TTCP.html for de¬ 

tails). This instrument allows the simultaneous observation 
of the variable and comparison stars, as well as of the sky. 
Therefore it is possible to obtain high precision photometry 
with fast sampling. The data reduction is also straitforward. 

2. The BAO measurements were obtained with the 0.85m tele¬ 
scope placed at the Beijing Astronomical Observatory. The 
telescope is equipped with a Three Channel Photometer 
(TCP) which is very similar to the TTCP quoted above. 
More details on BAO instruments can be found at 
www.bao.ac.cn. 

3. The OSN measurements were obtained with the 0.90 m 
telescope located at the Observatorio de Sierra Nevada in 
Granada, Spain. The telescope is equipped with a simul¬ 
taneous four-channel photometer (uvby Stromgren photo¬ 
electric photometer). More details on the instrument can be 
found at http://www.osn.iaa.es/fotometro^ng.html 

4. The SPM measurements were carried out with the 1.5 m 
telescope at the San Pedro Martir Observatory, Mexico. 
The 1.5 m telescope was equipped with a simulta¬ 
neous uvby Stromgren photoelectric photometer, which 
is the twin of the OSN instrument mentioned above. 
The SPM instrument (named “Danish”) is described at 
http;//haro.astrossp.unam.mx/Instruments/ 

5. The Serra la Nave data were obtained with the 0.91cm 
telescope of the Catania Astrophysical Observatory, 
Italy. The telescope is equipped with a Single Channel 
Photometer (SCP here after) which is described at 
http://w3c.ct.astro.it/sln/strumenti.html 

6. The APT measurements were obtained with the T6 
0.75 m Vienna Automatic Photoelectric Telescope 
(APT), situated at Washington Camp in Arizona, USA 
( [Strassmeier et al. 1997jBreger & Hiesberger, 1999^ . Also 
this telescope is equipped with a SCP. 

7. The SARA measurements were collected by using the 0.9 
m automated telescope at Kitt Peak in Arizona (USA). 
Images were collected using an Apogee AP7p camera with 
a back-illuminated SITe SIA 502AB 512x512 pixel CCD. 
The pixels are 24 microns square, corresponding to 0.73” 
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at the telescope focal plane scale. Sky flats, dark and bias 
exposures were taken every night. All data was calibrated 
and reduced using standard IRAF routines. 

8. The OGS observations were collected with the 1.0m 
OGS telescope placed at the Teide Observatory (Tenerife, 
Canary Islands, Spain). The telescope was equipped with a 
1024x1024 pixel CCD. The pixel scale was 0.32 for a to¬ 
tal field of view of 5.5 x5.5 . Typical integration times were 
10-15 sec, depending upon sky conditions. Sky flats, dark 
and bias exposures were taken every night. All data was 
calibrated and reduced using standard IRAF routines. 

9. The SOAO data were obtained with the 0.61m telescope 
placed at the Sobaeksan Optical Astronomy Observatory 
(SOAO) in Korea. The telescope was equipped with a SITe 
2K CCD camera. The field of view of a CCD image is 
about 20.5 x 20.5 arcmin^ on the f/13.5 Cassegrain focus 
of the telescope. Data reduction has been performed using 
the standard IRAF/CCDRED package. Instrumental mag¬ 
nitudes were obtained via the simple aperture photometric 
routine in the IRAF/APPHOT package. 



